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■ We use the model of ballistic deposition as a simple way to establish cooperation among the 

columns of a growing surface, the single individual of the same society. We show that cooperation 
, generates memory properties and at same time non-Poisson renewal events. The variable generating 

memory can be regarded as the velocity of a particle driven by a bath with the same time scale, 
and the variable generating renewal processes is the corresponding diffusional coordinate. 
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PACS numbers: 05.40.-a, 05.65.+b, 89.75.Da 



I. INTRODUCTION 



The Science of Complexity |1| is an emerging field of research, which is attracting an increasing number of investiga- 
tors. There is a general agreement that the deviation from the exponential behavior is a sign of complexity. However, 
on the origin of this deviation is there are different proposals. There are authors emphasizing the memory properties 
0,0,3 an d other 0, IE SB IS 01 stressing the non-Poisson renewal properties of complex systems. 

As an attractive example of non-Poisson renewal processes we refer to the field of single-molecule spectroscopy 
[TH and blinking quantum dots There is an increasing evidence that phenomena such as intermittent 
fluorescence yield histograms of the "on" and "off" sojourn time distribution that depart dramatically from the 
exponential condition. At the same time, a careful statistical analysis proves that these are renewal processes 0, ^| . 
Renewal theory ^(| describes, for example, successive replacements of light bulbs: when a bulb fails it is immediately 
replaced by a new one, which works with no memory whatsoever of the time duration of the earlier light bulb. Thus, 
t-H , renewal processes are characterized by the occurrence of events that reset to zero the system's memory. In which 
sense, therefore, a non-Poisson renewal process is a manifestation of organized behavior, which seems to imply 0,0,9 
*^ , long-time memory? 

■ It is tempting to express the non-exponential properties of a complex system in terms of ordinary Poisson processes 
with a fluctuating Poisson parameter 0,0,0. Using this perspective, it is possible to prove that complex 
_^ ■ networks emerge from fluctuating random graphs. However, it is worth remarking that slow fluctuations |2l| reduce the 
\Q ' occurrence of renewal events and generate memory properties compatible with the adoption of stationary correlation 

, functions. 

Another example of non-Poisson renewal process is given by the fields of seismic fluctuations l22l. where, however, 
there is no general agreement on the renewal character of the process. The authors of Ref. |22j| argue that the 
distribution of time distances between two consecutive main shocks corresponds to the properties of a non-Poisson 
renewal process. Is this conclusion reliable, or does it conflict with the notion of main shocks being complex processes, 
are |23|. to some extent, predicable? Another example, of greater interest for this paper, is given by the random growth 
of surfaces. The authors of a recent work [24j have used the concepts of non-Poisson renewal theory to describe and 
derive the same pro per ties that other authors 0,0,0] derive from the adoption of the Fractional Brownian Motion 
(FBM) perspective [28(. On one side, we have renewal, with jumps resetting to zero the system's memory, and on the 
• i-h , other we have fluctuations with infinite memory. Which is the best signature of complexity? 

' The purpose of this paper is to prove that both aspects are present in the processes of random growth of interfaces, 
which are a natural realization of of the Renewal Cooperation (RC) perspective proposed in the earlier work of Ref. 
29]. 

The outline of the paper is as follows. In Section [H] we shall give detail on the model used for the illustration of 
our approach to complexity, and we shall define the variable £(t) which is the signature of memory, and the variable 
y, whose origin regression is a non-Poisson renewal process. In Section IlIII we establish the renewal character of the 
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variable y and in Section Hvl we study the memory properties of the variable In Section [V] we shall point out the 
result of this paper. 

II. BALLISTIC MODEL 

We plan to shed light into the issues mentioned in Section H] with the model of Ballistic Deposition (BD) , which, 
is, in fact, a paradigmatic model for complexity. We shall limit our discussion to the (l+l)-dimcnsion case, and as 
a consequence, Fig.l is fully adequate to illustrate how the model works. At any time step, n = 1,2, •• • we select 
randomly one of the L columns, and we drop a particle on it. The particles fall down in sequence till to settle either 
at the bottom of the column or at the top of an earlier selected particle that by chance fell down in the same column. 
However, if the height of one of the two nearest neighbor columns is higher than the selected column, the particle sticks 
to the side of the highest particle of this neighbor column. There is also periodic boundary conditions to decrease the 
effect of limited size of surface. 

Actually, this side sticking action corresponds to a transverse transport of information, through which the column 
under study is informed about the height of the surrounding columns. Thus, examining the time evolution of a single 
column is equivalent to studying the behavior of a single individual and to assessing to what an extent it reflects the 
properties of the whole society. We plan to prove that the cooperation among the different columns generates memory 
and this memory generates renewal effects. 

It is well known |30| that the time distance between the arrival of one particle in this column and the next has the 
Poisson time distribution 

V>(r) = Aexp(-Ar), (1) 

where 

A =± (2) 

Of course, when one particle arrives, the height of the column increase by a quantity that can be also much greater 
than 1, thank to the side-sticking effect. We denote by the height increase of the column, and this quantity is 
when no particle arrives, and a number equal to 1 or larger when a particle arrives. 
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FIG. 1: Model of ballistic deposition. The particle B settles at the top of an earlier particle of the ame coloum, given the fact 
that there are no particles at a higher level in the two nearest neighbor columns. The particle A sticks to the right side of the 
left nearest neighbor colomn rather than at the top of a particle in the same colomn at a lower level. 
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The height of the column at a give time t is given by 



n=0 



with £o = 0. To produce a renewal process we use the variable y(t) defined by 

y{t) = h(t)- < h{t) >, 



(3) 



(4) 



where < h(t) > denotes the average over all the sample columns. In the earlier work of Ref. [2J| it was argued 
that the regressions of the variable y to the origin y = is a renewal process. We shall prove this property with a 
compelling numerical experiment. 

Let us now define the variable that is the signature of the memory properties generated by the BD model. This is 
the variable defined by 



(5) 



where the symbol £ denotes a time average. We generate a very long, but finite, sequence £(t), and we define the time 
average 



(6) 




10000 



FIG. 2: Probability distribution of jumps, f, (+), fitting to dashed line 0.22 exp(— 0.35 t); and waiting times of jumps, (x), 
fitting to solid line 0.001 ezp(-0.001 t). L = 1000 



Fig. 2 illustrates the distribution density tp(r) and and prove that the latter distribution, as well as the 

former, is an exponential, 



^>(£) = aexp(-6£). 



(7) 



On the basis of these results one would be tempted to conclude that the single column process is Poisson. We 
shall prove that it is not so because the fluctuation £(t) has memory, in spite of its exponential distribution. This 
memory is a consequence of the process of transverse transport of information, and it can be considered a signature 
of cooperation. 
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III. AGING 



The authors of Ref. [24| studied the variable y, whose variance w(L, t), in accordance with the literature, was shown 
to obey the anomalous prescription 



w(t) oc t , 



(8) 



with (3 < 0.5. Actually, for L — > oo the power index (3 is expected to fit the Kardar-Parisi-Zhang (KPZ) prediction [31j 
(3 = 1 /3. To speed up the numerical calculation we limit our calculations to the case of L = 1000, and consequently to 
the value (3 — 0.28. The authors of Ref. |2J] argued that the y-process can be reproduced with remarkable accuracy 
by adopting the subordination approach |5l Irl vh. I§l. ITof . It is important to stress that the variance w(t) is a Gibbs 
property, obtained by making a statistical average on all the columns of the sample. The single trajectory y(t) is 
derived from the stochastic trajectory of y(n), which is driven in the natural time scale by the ordinary Langevin 
equation 



dn 



y = -jy(ri) + f(n), 



(9) 



and assuming that the transition from the natural time scale n to the i-time scale is realized with the prescription 
t(n + 1) — t(n) = r(n), where r(n) is randomly derived from a distribution density with power law /i < 2. This is 
an approach well distinct from the ordinary memory approach 0, 0> 3 , which is based on a generalized Langevin 
equation, whose memory kernel is a conventional equilibrium correlation function. It is also important to point out 
that there exists a connection between /i and /3 stemming from It is worth remarking that all this is also a clear 
evidence of the fact that the relationship 



(3 = 2- ^ 



(10) 



which was derived by the authors of Ref. [24| from the renewal theory, rather than by using the FBM theory, as done 
by earlier investigators 1321 . 

The authors of Ref. [24| established a very good agreement between theory and numerical simulation by adopting 
this subordination perspective. In this section we confirm the arguments of Ref.|24| by means of an aging experiment, 
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FIG. 3: Probability distribution of waiting times of y(t) are shown with (+) when there is no aging. The straight solid line 
is 20 x" 1 ' 72 . Probability distribution of aged waiting times are dashed curve for experimental and solid curve for theoretical 
calculation. 
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recently proposed to prove the renewal character of BQD processes 15]. According to this procedure, we study first 
the recursions to the origin of the fluctuation y(t), which are described by Fig. 3. We see that y(t) remains in the 
positive portion of the y-axis for a given amount of time, then y moves to the negative portion of this axis, remains 
there for a while, till to re-cross the origin y = 0, and so on. We evaluate the histogram and we find an inverse power 
law with /i = 1.72. According to the prescription of Ref. ^^] we evaluate also the aged waiting time distribution 
ipi ex p) (t, i ). This is done as follows. The re-crossing times generates the time series {U}. For each time U we record 
the first time of the sequence at a distance from ti equal to or larger than ti + t a . This time will be t/-, with k > i. 
We record the time distance r(t,,t a ) = tk — (ti +t a ). We repeat the procedure for all the times of the sequence {ti}, 
and use the sequence of these recorded time distances to generate the distribution density ip^ exp ^ (r, t a ). We have also 
to evaluate 

^(r,t a ) = %^, (11) 
is. (t a ) 

where K(t a ) is a suitable normalization constant and ij)(t) is the experimental waiting time distribution corresponding 
to t a = 0. The accordance between ^ theor ^ (r,t a ) and ip( expS > (r,t a ) is judged Jl5j to be the numerical evidence of the 
renewal nature of the process under study. 

Fig. 3 confirms the renewal nature of the y-process, thereby providing further support for the theoretical approach 
to the stochastic growth of a single column proposed in Ref. [24J and for the adoption of the renewal theory to derive 
the crucial relation of Eq. (|10f> . Using Eq. (|10fl and the numerical results of Fig. 3 for fj,, we obtain f3 = 0.28. This 
same power index will be derived from the memory properties discussed in Section llVl 



IV. MEMORY 



In this section we adopt a single-trajectory perspective, rather than the Gibbs approach of Section flHI How can we 
reveal the cooperative nature of the process if we limit ourselves to studying the time evolution of only one column? 
Apparently, the statistical analysis that led us to Eqs. JIJ and Q would suggest that the anomalous growth property 
of the interface as a whole is annihilated by the observation of a single column. 

To assess the memory properties created by cooperation, we proceed as follows. We create the diffusing variable 

x(t) = /%"(*')• (12) 
Jo 

Then, using the proposal made by Stanley and co-workers years ago [33, 031 we convert the single diffusing trajectory 
of Eq. (|12|l into many diffusing trajectories. This is done as follows. We consider a window of size I, and we move 
it along the sequence. We record the space positions of the random walker at times s and s + I. This makes it 
possible for us to define the s-th random walker that in the time I, moves from the origin x — to the position 
x(s, I) = x(s + I) — x(l). Thus we have a set of walkers, each walker denoted by the label s, moving for a time / and 
covering the distance x(s,l). We imagine all these walkers to depart from the origin x = at the same time. This 
means that the set of values x(s, I) determines a distribution, whose variance w(l) is given by 

w(l) = ((x( S ,iy<x( S J)>) 2 y /2 , (13) 

the averages < ... > being made on all the labels s. It is evident that if the single column reflects the complex behavior 
of the whole surface, see Eq. (JSJ, then 

lim log(w(/)) =/31og(Z). (14) 

l — >oc 

Fig. 4 shows that this prediction is fulfilled and that, as expected, (3 — 0.28. 

Under the stationary assumption, it is straightforward to prove (see, for instance, Ref. |35j) that 

j t <x 2 {t)>=<e> J q ^df, (15) 

where denotes the equilibrium correlation function of the fluctuation £(i). In the time asymptotic limit, which 

is the intermediate regime illustrated by both Fig. 3 and Fig. 4, with no upper bound time limit, determined by the 
third regime, we get 



w (t) 2 =< x 2 {t) >ozt 2(} , 



(16) 




t 

FIG. 4: Variance of x is shown with dashed line. The solid line is 0.256 a; 0,288 



namely, we reiterate that the diffusion process has anomalous scaling defined by /3, with (3 < 1/2. By differentiating 
Eq. H15fl with respect to time we establish a connection with the time asymptotic properties of ^^(t). Due to the fact 
that P < 1/2 we conclude that the correlation function must have a negative tail, namely that 

_ , . constant , 
ton* f (*) = — (17) 

with 



ip = 2 - 2/3. (18) 

This property establishes a closer connection with the earlier work of Ref. 12911 . where the variable responsible for 
memory yields a correlation function with a negative tail. In the case of Ref. |29| the correlation function is known 
theoretically, so that it is possible to move from the correlation function to the variance time evolution using Eq. JTSj. 
In the present case, we do not have available any analytical approach to the equilibrium correlation function of £(t). 
Its numerical derivation, as pointed out by Eq. (|15fl . would be equivalent to differentiating twice the variance w(t), 
which is numerically a source of big errors. In, fact, the numerical approach to the equilibrium correlation function, 
not shown here, yields a negative tail, but the assessment of the correct power requires rich statistics and excessive 
computational time. 

Thus, to do double check our conclusion, that the variable £(t) is the source of memory, with no renewal events, 
we prefer to accurately examine the physically plausible conjecture made by the authors of Ref. 36] that the non- 
Poisson renewal events may be hidden among the pseudo-Poissonian fluctuations of the variable According to 

this conjecture, either the fluctuations of above a given threshold R or the time distances between two consecutive 
fluctuations of exceeding the threshold R, may deviate from the exponential distribution. Fig. 5 and Fig. 6 
illustrate these two conditions, and prove that in both cases the asymptotic time behavior is exponential. 

We therefore conclude that there is no evidence of renewal events in the dynamics of £(t). The stationary assumption 
made to deal with the correlation function of £(t) does not conflict with the non-stationary nature of y. This is so 
because £(f) and y(t) are the signatures of cooperation and renewal, respectively. Cooperation generates a slow 
deca ying correlation function, but the diffusion process generated by these fluctuations is renewal, in accordance with 
Ref. Hf. 
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V. CONCLUDING REMARKS 



This paper shows how to reconcile a perspective based on memory with one based on renewal non-Poisson processes. 
To a hrst sight, these two visions can be perceived as bein g in compatible. For instance, the adoption of the FBM 
theory adopted by some authors to account for persistency l37l [38l l3^| seems to conflict with the renewal approach. 
On the contrary, in accordance with the recent result of Ref. [23 this paper shows that the cooperative nature of the 
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BD model generate memory and, consequently, anomalous diffusion. However, the origin re-crossings of the variable 
y are renewal, in the same way as the diffusion variable generated by the long-memory fluctuations used to generate 
dynamically FBM exhibit a renewal character . 

We acknowledge financial support from Welch foundation through Grant N. 70525 and from ARO through Grant 
No. W911NF-05-1-0059. 



[1] R. Badii and A. Politi, Complexity -.Hierarchical Structure and Scaling in Physics (Cambridge University Press, Cambridge, 
UK, 1997). 

[2] U. Balucani, M. H. Lee, V. Tognetti, Phys. Rep. 373, 409 (2003). 
[3] S.S. Melnyk, O.V. Usatenko, V.A. Yampl'skii, Physica A 361, 405 (2006). 
[4] W. Min, G. Luo, B. J. Cherayil, S.C. Kou, and X. S. Xie, Phys. Rev. Lett. 94, 198302 (2005). 
[5] LM. Sokolov, Phys. Rev. E 63, 011104 (2000). 
[6] I. M. Sokolov, Phys. Rev. E 63, 056111 (2001). 
[7] I. M. Sokolov, Phys. Rev. E 68, 041101 (2002). 
[8] E. Barkai and R.J. Silbey, J. Phys. Chem. B 104, 3866 (2000). 
[9] R. Metzler and J. Klafter, J. Phys. Chem. B 104, 3851 (2000). 
[10] A. Budini, Phys. Rev. A 69, 042107 (2004). 

[11] E. Barkai, Y. Jung, and R. Silbey, Annu. Rev. Phys. Chem. 55, 457 (2004). 

[12] M. Nirmal, B. O. Dabbousi, M. G. Bawendi, J. J. Macklin, J. K. Trautman, T.D. Harris, L.E. Bras, Nature 383, 802 
(2006). 

[13] M. Kuno, D.P. Fromm, H.F. Hamann, A. Gallagher, D.J. Nesbitt, J. Chem. Phys. 115, 1028 (2001). 

[14] X. Brokmann, J. -P. Hermier, G. Messin, P. Desbiolles, J.-P. Bouchaud, M. Dahan, Phys. Rev. Lett. 90, 120601 (2003). 

[15] S. Bianco, P. Grigolini, P. Paradisi, J. Chem. Phys. 123, 174704 (2005). 

[16] DR. Cox, Renewal Theory (Wiley, New York, 1962). 

[17] G. Wilk and Z. Wlodarczyk, Phys. Rev. Lett. 84, 2770 (2000). 

[18] C. Beck, Phys. Rev. Lett. 87, 180601 (2001). 

[19] C. Beck and E.G.D. Cohen, Physica A 332, 267 (2003). 

[20] S. Abe and S. Thurner, Phys. Rev. E 72, 036102 (2005). 

[21] P. Allegrini, F. Barbi, P. Grigolini, and P. Paradisi, Phys. Rev. E 73, 046136 (2006). 

[22] M. S. Mega, P. Allegrini, P. Grigolini, V. Latora, L. Palatella, A. Rapisarda, and S. Vinciguerra, Phys. Rev. Lett. 90, 
188501(2003). 

[23] X. Yang, S. Du, and J. Ma, Phys. Rev. Lett. 92, 228501 (2001). 

[24] R. Failla, M. Ignaccolo, P. Grigolini and A. Schwettmann, Phys. Rev E 70 R, 010101 (2004). 

[25] J. Krug, H. Kallabis, S.N. Majumdar, S. J. Cornell, A. J. Bray and C. Sire, Phys. Rev. E 56, 2702 (1997). 

[26] J. Merikoski, J. Maunuksela, M. Myllys, J. Timonen, and M. J. Alava, Phys. Rev. Lett. 90, 024501 (2003). 

[27] S.N. Majumdar and D. Das, Phys. Rev. E 71, 036129 (2005). 

[28] B. B. Mandelbrot, The Fractal Geometry of Nature (Freeman, New York, 1998). 

[29] R. Cakir, P. Grigolini, A. Krockin, Phys. Rev. E, in press (2006). 

[30] H.E. Stanley and A. L. Barabasi, Fractal Concepts in Surface Growth, Cambridge University Press, Cambridge (1995). 
[31] M. Kardar, G. Parisi, and Y.-C. Zhang, Phys. Rev. Lett. 56, 889 (1986). 
[32] M. Ding and W.Yang, Phys. Rev. E 52, 207 (1995). 

[33] C.K. Peng, S.V. Buldyrev, A.L. Godberger, S. Havlin, P. Sciortino, M. Simons, and H.E. Stanley, Nature (London) 356, 
168 (1992). 

[34] Z. Chen, P.Ch. Ivanov, K. Hu, H. E. Stanley, Phys. Rev. E 65, 041107 (2002). 
[35] G. Trefan, E. Floriani, B. J. West and P. Grigolini, Phys. Rev. E 50, 2564 (1994). 
[36] R. Failla and P. Grigolini, Fluctuation and Noise Letters, 5, L175 (2005). 
[37] S. N. Majumdar and D. Das, Phys. Rev. E 71, 036129 (2005). 

[38] M. Constantin, C. Dasgupta, P. P. Chatraphorn, S. N. Majumdar, and S. Das Sarma, Phys. Rev. E 69, 061608 (2004). 
[39] J. -F. Muzy, E. Bacry, Phys. Rev. E 66, 056121 (2002). 



